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Abstract. The observational constraints on the primordial power spectrum have 
tightened considerably with the release of the first year analysis of the WMAP 
observations, especially when combined with the results from other CMB experiments 
and galaxy redshift surveys. These observations allow us to constrain the physics of 
cosmological inflation: 

(i) The data show that the Hubble distance is almost constant during inflation. 
While observable modes cross the Hubble scale, it changes by less than 3% during 
one e-folding: dn < 0.032 at 2f7. The distance scale of inflation itself remains 
poorly constrained: 1.2 x 10^^^ cm < dn < 1 cm. 

(ii) We present a new classification of single-field inflationary scenarios (including 
scenarios beyond slow-roll inflation), based on physical criteria, namely the 
behaviour of the kinetic and total energy densities of the inflaton field. The 
current data show no preference for any of the scenarios. 

(iii) For the first time the slow- roll assumption could be dropped from the data 
analysis and replaced by the more general assumption that the Hubble scale 
is (almost) constant during the observable part of inflation. We present simple 
analytic expressions for the scalar and tensor power spectra for this very general 
class of inflation models and test their accuracy. 
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1. Introduction and results 

With the release of the analysis of the first year WMAP data ^1 E] and the first 3d 
SDSS power spectrum jS], we got splendid confirmation of the long-standing expectation 
that the primordial power spectrum of density fluctuations is almost scale- invariant, 
as anticipated already by Harrison and Zel'dovich in the 1970s. The simplest and 
most elegant mechanism known to produce primordial spectra with that property is 
cosmological inflation. Its beauty lies in the fact that we do not need to invent 
untested physical principles: quantum mechanics and general relativity in four space- 
time dimensions are good enough. What is needed is a cosmic substratum that gives 
rise to an epoch of accelerated expansion of the very early Universe. A positive energy 
density of the vacuum is the simplest example. Trying to figure out the underlying 
mechanism leads to the pressing question: What is the scale of cosmological inflation? 
It is probably the most poorly constrained energy scale in cosmology; we could argue 
that cosmological inflation happens somewhere between the electroweak scale (requiring 
a mechanism to create baryons after the end of inflation) and the highest possible energy 
scale, which might be the scale of quantum gravity. This implies an uncertainty of some 
16 orders of magnitude. 

Although we have no clue on how the world looks like at scales beyond the standard 
model of particle physics, we can make some generic predictions about the fluctuations 
in the energy density and space-time structure that are seeded by quantum fluctuations 
during cosmological inflation jUEj- 

We assume that the first prediction of cosmological inflation, spatial flatness of 
the Universe, is well established and take Q = 1 (see [0112] for a detailed discussion). 
Further, we assume that the observed perturbations are of isentropic nature, which is 
another prediction of cosmological inflation, unless there is some matter component 
in the Universe that never ever coupled to the radiation fluid. The most convincing 
evidence for the dominance of isentropic fluctuations comes from the recent detection 
of polarization in the cosmic microwave background (CMB) [ZllHl, which requires the 
presence of a quadrupole anisotropy at the moment of photon decoupling and therefore 
shows that large-scale fluctuations existed back then. In the following we focus on the 
power spectra of scalar and tensor fluctuations, which contain all the information if, 
again as predicted by inflation, the fluctuations are Gaussian (which is consistent with 
the data [Hj). 

Two important inputs are required to predict the inflationary power spectra: the 
Hubble rate H during cosmological inflation as a function of the scale factor a (or, for 
later convenience, as a function of the logarithm of the scale factor, denoted by N) and 
the number of dynamical degrees of freedom that drive inflation. In the case of de Sitter 
space-time this number is zero, but in such a model inflation does not end; the vacuum 
energy thus must be dynamical. In the simplest scenarios with an exit from inflation 
a single scalar field if, usually called inflaton, is responsible for the acceleration of the 
Universe. We will restrict the following considerations to that case, with equations of 
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motion, 



(p + m^ + V'{^) = Q, (1) 
1 
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(2) 



Here, V{ip) is the potential energy density of the inflaton field, the dot stands for a 
derivative with respect to cosmic time t, and the prime denotes a derivative with respect 
to the infiaton field. We use Planck units, SttG = c = h = 1. 

In principle, observations should allow the reconstruction of H{N), but since only 
a finite interval of wave numbers k is accessible to observations, we can only probe a 
finite and rather short interval AN. An efficient encoding of H{N) is by the so-called 
horizon-fiow functions, evaluated at a pivot point adapted to the experimental set- 
up. The horizon-fiow functions are a generalization of the slow-roll parameters ^Ul and 
are defined recursively as the logarithmic derivatives of the Hubble scale with respect 
to the number of e-foldings N 

Wi-^^ Vm>0, eo = ^, (3) 

where iVj denotes an arbitrary 'initial' moment. The necessary condition for infiation 
to take place (d > 0) becomes ei < 1, and we assume here that the weak energy 
condition and null energy condition hold true, i.e. ei > 0. We define the graceful exit 
from infiation as the moment when ei crosses unity. Specifying the set {em(A^*)} is 
equivalent to specifying H{N). A truncation of the set {em(A^*)} corresponds to an 
incomplete knowledge of the evolution of the Hubble rate. 

In single field slow-roll models of infiation, all the horizon-fiow functions are 
typically small and, hence, equivalent to the slow-roll parameters. In terms of the 
infiaton potential V{ip) and its derivatives with respect to the infiaton field, the first 
two fiow functions are given by fT^ : 




^2 




(4) 



CMB experiments (the most recent results come from CBI [1^1 116j. Archeops |14j . 
ACBAR Hg, VSA [milHI and WMAP P) and galaxy redshift surveys (especially 
2dF [iniEII! and SDSS |3]) have measured the amplitude A and the spectral index n of 
density fiuctuations and constrained the running^ of the spectral index dra/d In k, as well 
as the amount of gravitational waves, typically expressed as the tensor-to-scalar ratio 
r, that can contribute to the CMB signal. This set of observables {A, n,dn/ din k,r} 
could, in principle, provide a measurement of the infiationary parameters {H, ei, 62, €3}. 

In a number of works, various combinations of experiments and methods have been 
used to determine a subset of these parameters, see especially Refs. pill^l^l^l^l^ . 

I Two recent works on new data from CBI |TB1 and VSA claim to see evidence for a large (compared 
to expectations from cosmological inflation) running of the spectral index with negative sign, but point 
out a calibration issue that weakens the evidence. 
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Most interesting are the constraints in the €1-62 plane; see Fig. The common finding 
is that models with almost scale-invariant spectra provide acceptable fits and ei < 0.032 
(at 2a) |21]. The corresponding likelihood contours from that work are shown in Fig. ^ 
where the dashed line denotes scale-invariance, up to third order corrections in the 
horizon-flow functions. Thus, observations show for the flrst time that ei <^ 1. The 
limits on €2 are much less restrictive; the situation |e2| > ei covers a large part of the 
allowed parameter space. This is due to the so-called tensor degeneracy ^7\, i.e. one can 
compensate an increase of the tensor contribution by making the scalar spectral index 
bluer. Since there exist only upper limits on the contribution from gravitational waves, 
the scale of inflation cannot be flxed by present observations. However, the combination 
of the upper limit on ei and the measurement of the amplitude of scalar fluctuations 
also allows an upper limit on the energy scale of inflation V^^^ < 0.01 (= 2.7 x 10^^ GeV) 
j24j or a lower limit on the distance scale during inflation (in > 1-2 x 10^^^ cm. A lower 
(upper) limit on the energy (distance) scale can be obtained from the facts that the 
Universe contains baryons and that there is no known mechanism to produce baryons 
below the electroweak scale (100 GeV, respectively 1 cm) [2H]- Thus observations put 
the scale of inflation at least two orders of magnitude below the Planck scale [this limit 
will improve proportionally to the limits on ^/ii and will be tightened as the value of 
the optical depth r is pinned down more precisely, since Aexp(— 2r) is the observed 
quantity]. 

Based on the observational constraints on the flrst two horizon-flow functions, we do 
not restrict our considerations to slow-roll models in this paper. We relax the slow-roll 
conditions (e^ ^ 1 for all m) to < 1 for some m. In Section |21 we demonstrate that 
the condition ei <^ 1 is sufficient to guarantee that the Hubble scale is almost constant 
during the infiationary epoch of interest. 

The physical meaning of the first and second horizon-fiow functions is discussed 
in Section IHl in the framework of single-field models, and a new physical classification 
of infiationary scenarios is introduced, based on the behaviour of the kinetic and total 
energy densities of the infiaton. To our surprise this approach closely resembles the 
small field/large field/hybrid classification of Dodelson et al. [22] 5 but is not restricted 
to slow-roll models. 

In Section 4 we derive a procedure to self-consistently calculate higher-order 
corrections at any scale k, given an analytic expression of the power at the pivot scale 
fc* as a function of the horizon-fiow functions. A similar relation has been used in the 
literature to calculate the spectral index and its running so far, but, to our knowledge, 
it has never been justified rigorously. 

Starting from the observation that ei <C 1, we show that the slow- roll approximation 
is no longer needed in the analysis of the data and that it could be replaced by a 
more efficient approximation (constant-horizon approximation JI]). The latter is more 
efficient in the sense that we need to know less horizon-fiow functions in order to predict 
the power spectra to a given accuracy, as compared with the slow-roll approximation. 

We map the classification of models from the €1-62 plane to the space spanned 
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hidden-exit inflation 



toward-exit inflation 




Figure 1. Observational constraints, interpretation, and accuracy of various 
approximations in the £1-62 plane: the likelihood curves have been provided by 
S. M. Leach and enclose the la-, 2a- and Su-allowed regions. Models with a 
scale- invariant spectrum fall on the dashed line (HZ). An observation of 62 > would 
favour simple models of inflation, in which the evolution toward a graceful exit is seen 
(toward-exit inflation). If £2 < 0, a more contrived scenario must be realized, since 
in that case the mechanism of exit is hidden from observation (hidden-exit inflation). 
The region with £2 > is further split up by the (full straight) line £1 = £2/2, which 
distinguishes models in which the kinetic energy density decreases (above the line) or 
increases with time. The latter requires a mechanism to ensure that the kinetic energy 
density is very small initially. The solid arc centered at the origin encloses the region in 
which the slow-roll approximation at second order (sra2) provides a prediction of the 
amplitude at the pivot scale better than 0.3%. Below the hyperbolic lines the same 
accuracy is achieved for the constant-horizon approximation (cha2 at second order, 
cha3 at third order) and the growing horizon approximation (gha). The dotted lines 
represent the boundaries between cha3 and cha2, as well as between cha2 and gha. 



by the tilt of the spectrum, n — 1, and the ratio of tensor and scalar amplitudes r in 
Section 5. At leading order in the horizon-flow functions this is a one-to-one mapping; 
including higher orders introduces an ambiguity, due to the running of the spectral 
index. We thus argue that the €1-62 plane is more fundamental. For practical purposes, 
the r-tilt plane is nevertheless a useful tool. 

In Section 6 we give examples to demonstrate the accuracy of the new method and 
discuss the applicability of the new method in the light of the recent data. 

Our main results are summarized in Fig. ^ We overlay our new classification of 
models on the likelihood contours as obtained by Leach and Liddle P^. We also indicate 
which approximation to the power spectrum at the pivot scale is best suited for which 
region of parameter space. We assume an accuracy goal of 0.3% at the pivot, which puts 
the theoretical errors safely away from the systematic and statistical errors of WMAP 
and Planck. 
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2. Hubble horizon flow during inflation 

The Hubble horizon, dn = c/H (recall that c = 1), is roughly the size of the region where 
causal processes can take place during one Hubble time 1/H. An inflationary epoch is 
characterized by a decreasing comoving Hubble horizon, dn/a. Constant vacuum energy 
is the simplest (but unphysical, since inflation would continue forever) kind of matter 
leading to an inflationary universe. In this scenario, the Hubble scale is constant. More 
generally, during inflation, dn is expected to vary slowly during a given number of e- 
foldings N: dN/dt = H. The behaviour of the Hubble horizon can be parameterized 
as 



where the coefficients in this expansion are expressed in terms of the horizon-flow 
functions at the pivot point. One can see that, if ei ^ 1 and < 1 for any 
m > 1, then ~ constant. For such cases, the smaller ei is, the larger the other 
horizon-flow functions can be, since ei is the leading factor of all coefficients of the 
higher-order terms in the series. We can use this observation as the starting point for 
what we call the constant-horizon approximation (ei ^ 1 and < 1 for all m > 1). 
As is seen in Fig. [T] and discussed in the introduction, the data show that c^h = ci <^ 1. 

3. Classification of single scalar field inflationary models 

The number of proposed models of cosmological inflation is very large and keeps growing. 
Current observations allow us to constrain the number of successful models, and future 
observations are expected to play a stronger role in that direction. For a more eflective 
use of the data analysis, it is useful to group the inflationary scenarios according to 
some commonly applicable criteria, i.e. independent of parameters such as masses or 
coupling constants, which are strongly mo del- dependent. For slow-roll inflation, such a 
classiflcation was suggested by Dodelson, Kinney and Kolb [211, based on the shape of 
the inflaton potential (and expressed as conditions on the slow-roll parameters). 

Here we present a new scheme in terms of the horizon-flow functions, which closely 
resembles the Dodelson et al. classiflcation, but avoids misleading terminology, allows 
the inclusion of models beyond the slow-roll class, and is based on physical criteria on 
the kinetic and total energy densities. 

The idea is to investigate how the kinetic and potential energy of the inflaton fleld 
change with time, both absolute and relative to each other. It is thus useful to rewrite 
the flrst two horizon-flow functions as 
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where we used Eqs. (H)) and 0. Using Eqs. (jT)) and ((Tj), the time variation of the 
potential energy is found to be 



Since Hip"^ > 0, the potential energy density can never increase for physically meaningful 
values of ei and €2- Thus, we focus on the behaviour of the kinetic energy density. 

According to Eq. (0), ei/3 > measures the ratio of kinetic energy density to total 
energy density. Using definition , we can ask how this ratio changes with time during 
inflation: 

^^=i/^6,. (9) 
dt 3 3 ^ ' 

Since H and ei are positive, €2 = marks a borderline between two physically different 

cases: increasing (e2 > 0) and decreasing (e2 < 0) kinetic energy density with respect 

to the total energy density of the inflaton. 

For a more complete understanding of the inflationary dynamics, we must also look 

at the absolute time variation of the kinetic energy density d(v9^/2)/d)f:, which, using 

Eqs. dTJ and (jH)), can be written as 



For ei > we have 7^ and thus €2 = 2ei is another borderline between two different 
physical scenarios in which kinetic energy density grows (e2 > 2ei) or falls (e2 < 2ei) 
with time. We thus arrive at the following classification: 

• €2 > 0: the kinetic energy density is increasing with respect to the total energy 
density. This is a necessary condition for evolving toward a graceful exit of inflation. 
Thus, such models could be called toward- exit models, since it can be argued that 
the approach to a graceful exit of inflation is observed in that case (although a more 
complicated interpretation remains possible). 

There are two subcases, corresponding to whether or not the kinetic energy density 
grows with time. It would be expected that this provides a criterion to discriminate 
between the false vacuum models and chaotic models in the slow-roll phase, since 
false vacuum models start out with vanishing kinetic energy density, whereas chaotic 
models always have a non- vanishing kinetic energy density. 

— 62 > 2ei: kinetic energy density grows with time. An example of this 
are false vacuum models arising, e.g. in superstrings models, with potential 

V = Vo — rn^ip^ /2 [SUj. Here, in the slow-roll phase, €2 ~ 4(Vo/m^(y9^)ei, and 
the condition is met as long as Vq > m?ip'^/2. False vacuum models need 
to provide a mechanism that gives rise to small kinetic energy density of the 
inflaton before observable modes cross the Hubble scale during inflation. 

— €2 < 2ei: kinetic energy density is decreasing with time. Monomial potentials 

V = Xif^/n with chaotic initial conditions give rise to models of that kind. 
During their slow-roll phase €2 ~ (4/'ri)ei, and the criterion is met for n > 2. 




(8) 




(10) 



8 



Here no extra mechanism is needed to tune the kinetic energy density initially, 
and one could argue that these are the simplest models. 
— €2 = 2ei: kinetic energy density is constant. For slow-roll models, this is 
realised for chaotic inflation with a quadratic potential (n = 2) pT] . 
Note that this criterion is different from the one given by Dodelson et al. [SH] to 
distinguish between small- and large-fleld models. In our notation and in the 
slow-roll approximation their borderline is €2 — 4ei. This is the case of a linear 
potential (n = 1). Naively one would think that these models should fall into 
the same class as chaotic inflationary models, but actually they do not give 
rise to a successful scenario, since for the linear potential V = Vo(l + v?/v'o) we 
flnd during slow- roll ei ^ €2/4:^ l/(2(y9Q), which is constant and thus inflation 
never ends. 

• €2 < 0: kinetic energy density is decreasing absolutely and relatively. Here, in 
order to reach a graceful exit, €2 has to change sign at some point. Thus, negative 
values of €2 must correspond to models in which inflation still has to go through 
a transition to either another stage of inflation or to directly to stop it by some 
unknown mechanism. The actual mechanism of exit is out of sight of observations, 
so we could call these models hidden-exit models. 

An example is the hybrid model, with an effective potential V = Vq + m^ip^ /2 in 
the slow-roll phase |H2]- In that case 62 ~ — 4(Vo/m^(y9^)ei. In order to end inflation 
Vo must depend on yet another fleld which flnally drives Vq to zero. 

• 62 = 0: the ratio of kinetic to total energy density is constant. If additionally 
all em = 0, then this is de Sitter space-time. If only ei 7^ 0, then this is power- 
law inflation jSS]. There also exists a number of models where 62 asymptotically 
converges to zero during inflation, so that they are observationally indistinguishable 
from power-law inflation (see [S^ for a discussion and further references). In all 
these cases there is no exit from inflation but the models are not ruled out by 
current cosmological data. 

Provided the slow-roll approximation holds, €2 = coincides with the limit between 
hybrid models and large-fleld models in the Dodelson et al. classification. 

Summarizing, our three classes are i) hidden-exit inflation (62 < 0), ii) toward-exit 
inflation with general ("chaotic") initial conditions (0 < £2 < 2ei), and iii) toward-exit 
inflation with special (e.g. false vacuum) initial conditions (0 < 2ei < £2). 

At this point it is convenient to note that the here introduced by us classiflcation 
involves only exact expressions. We confront this new classiflcation with the 
observational constraints in Fig. ^ Although the biggest piece of allowed parameter 
space falls into the hidden-exit inflation class, this cannot be seen as a preference of 
the data for this scenario, since the shape of the likelihood curves is due to the tensor 
degeneracy. At present all scenarios are consistent with the data. A better determination 
of the spectral index (lines parallel to the dashed line correspond roughly to flxed values 
of the spectral index) could in principle rule out some of the possibilities; there is thus 
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hope to learn more about inflation, well before we can expect to get a hand on the 
tensor contribution via observations of the B-polarization pattern of the CMB. 

4. Power spectra 

The accurate prediction of inflationary perturbations have been of concern since 
Mukhanov and Chibisov |1] realized that density and space-time perturbations during 
inflation could be the seeds for large-scale structure formation. The calculation of 
inflationary power spectra requires the solution of the mode equations for scalar and 
tensor fluctuations. The assumption that scalar and tensor perturbations are quantum 
fluctuations of the vacuum originally fixes the power spectra uniquely. The mode 
equations are of the same type as the Schroedinger equation, which happens to be 
difficult to solve, even for the simple models (see jSSl EHl IHZl IHH HH 1121 ^\ for details 
and references). 

It proved useful to expand the primordial power spectrum in a Taylor series around 
a pivot scale = {aH){N^:). Within a time interval At ~ l/H, the modes in the 
logarithmic frequency interval Aln/c = \n{eaH) — ln{aH) = 1 cross the Hubble scale§. 
It is thus natural to expand in terms of ln{k/k^,): 

V{k)=V{k.)j:^\n-{^), (11) 

n>0 "'■ 

where the coefficients are defined by 

d" V{k) 
dhFkViK] 

These coefficients depend on the horizon-flow parameters in such a way that they are 
regular in the limit — ^ 0; in fact oq — 1 and a„ — >■ for all n > 0. Now it becomes 
obvious why we have separated a factor that represents the leading-order prediction for 
the amplitude at the pivot scale 

-Psik.) ^ ^ (13) 
for scalar perturbations, and 

-PAh) ^ ^ (14) 

for tensor perturbations. In the notation of the WMAP team, Vs{k^) = 10~^y4 and 
VT{k.)/Vs{k.)=r. 

4-1. Independence from the pivot point 

The physical power spectrum must not depend on the arbitrary choice of a pivot scale 
k^, in the above expansion: 

§ The 'crossing' of the Hubble horizon is defined here to take place when k ^ aH. 



(12) 

k=k. 
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Evaluation of this expression provides us with the non-trivial relation 



n>0 ^- \^*/ 



dlnP 
d In fc, 



■an + 



din fc* 



-E 

n>l 



In 



n-l 



n 



- =0. (16) 



A comparison of the coefficients of In^lk/k^) finally leads to a recursion relation for the 
higher coefficients of the Taylor series: 



1-ei 



d\nV 



Vn > 



17) 



dN, din AT, 

where in the last step we used the identity 

dlnfc, = (l-ei)dAr,. (18) 

For scalars dlnPs/dlnA^* = — 2ei — £2, whereas for tensors dlnPT/dlniV* = — 2ei. Let 
us note that the recursion relation (jl7|) is an exact result: no approximation has been 
made, apart from the assumption that the linear perturbation analysis is justified. 



4.2. Approximation schemes 

Up to date, three approximation schemes have been proposed to allow for a high- 
precision calculation of the spectra amplitudes. Using recursion (|T7jl . all of these 
approximations yield expressions for the coefficients a„ as expansions in terms of the 
horizon-flow functions. Keeping terms up to order e'' (which stands here for any 
monomial of e^'s at order q) in Oq allows us to calculate all terms up to order e'^^"' 
in a„. We will denote the order of a given approximation by the highest order in the 
seed Oq. Besides the order q, a second choice that must be made prior to data analysis 
is how many coefficients a„ should be included in the analysis. The common practice 
is that only the terms Oq and ai are taken into account, 02 being included when the 
running of the spectral index is included. 



4-2.1. Slow-roll approximation. Assuming all horizon-flow functions to be small 
(without assuming any hierarchy among them), the equation of modes can be solved 
by an iterative method using Green's functions [3Hj. To second order in the slow-roll 
parameters, here expressed as horizon-flow parameters, the seed for recursion (fT7j) for 
the scalar spectrum reads 

aso = l-2(C+l)ei-Ce2 

+ (c'^-C + l!^-7),e.+ (44)..e3. (19) 

where C = 7E + ln2 — 2?« —0.7296, while for the tensor spectrum we have fT^ . 
aTo = l-2(C + l)ei 

+ (2C^ + 2C + ^-5]el+(-C^ -2C+^-2]eie2. (20) 
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In Fig. n we estimated the region of parameter space in which the prediction of the 
pivot amphtude A{k^) is better than 0.3% for the slow-roll approximation at second 
order. Such a high precision is actually needed to ensure that the power spectrum can 
be predicted with an accuracy better than a few per cent over at least three or four 
decades in wave number. 

4-2.2. Constant-horizon assumption. In a number of inflationary models the time 
derivative of the Hubble distance is tiny (see Section IHl for examples). For this kind 
of models, during a certain number of e- foldings, ci ^ 1. However, as we have seen in 
Section 121 this does not necessarily mean that all other have to be small as well. We 
thus worked out the constant-horizon approximation JT] at order q for the situation 
leal > max(|eie2|, |e2e3|), which means that, for oq we are allowed to include the following 
monomials in the primordial spectra: 1, ei, e2, . . . , With this approximation, the 
g = 3 expression for the seed of the scalar spectrum is 



Also here we estimate the region of 0.3% accuracy of A{k^) and indicate it in Fig. ^ 
For more details on how to do this we refer the reader to the work of |llj . 

4-2.3. Growing-horizon assumption. This approximation is valid for cases with |em| < 
ei for m > 1. The trivial example here is power-law inflation [S3], where only ei 7^ 0. 
This approximation is also valid in inflationary scenarios, where power-law inflation 
dynamics is a past or future attractor, and ei is sufficiently large for higher-order 
corrections to make sense. An expression for oq can be obtained by keeping all terms in 
ei up to order q, where q is the maximal integer for which el > max(|eie2|, |e2e3|) holds 
true. We defined the (linearly) growing-horizon approximation to order q to include 
the following terms: 1, ei, . . . , e^, £2. To third order the corresponding expression for the 
scalar spectrum is 




(21) 



where ({3) ~ 1.2021, and to any order for the tensor spectrum, 
aTo = l-2(C + l)ei. 



(22) 




(23) 



and for the tensor spectrum, 

axo = 1 - 2(C + l)ei 
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- (^ + C7r^-12C + ^C(3)-f)6?. (24) 
5. Classification of inflation models in the r— tilt plane 

In this section we indicate how to make contact with the variables that are frequently 
used in cosmological parameter estimation. However useful these variables have been 
assumed to be so far, as we shall see here, the analysis using them becomes much 
more involved that while using the horizon-flow functions. Complications arise because 
there are not exact model independent expressions for (n — 1) and r; they are given as 
expansions in terms of the horizon-flow functions. 

Using one or the other set of variables is just a matter of choosing the working 
parametrization for the primordial spectra. Very often a power-law shape is assumed 
for the scalar power spectrum; 

^s(^) = AiK) (^A^ , (25) 

where n is called the spectral index and n — 1 the tilt of the spectrum. A comparison 
with expression (jlll) reveals that 

A{h) = 'P{K)aso, (26) 

(n-l){K) = ^. (27) 
o-so 

For the slow-roll and the constant-horizon approximation we flnd at second order in the 
horizon-flow functions 

{n - = -2ei - 62 - 2el - (2C + 3)eie2 - Ceaeg. (28) 

A difference between the two approximations shows up at the third order. Here we 
restrict the discussion to the second order expressions. 

Tensor contributions are typically introduced via the tensor to scalar ratio 

r = ^4rV = 16e,(l + Ce2) (29) 

at second order in the slow-roll and the constant-horizon approximation. 

A closer inspection of the above expressions shows that there is a simple one-to-one 
map (ei, €2) ^ {r,n — 1) at the leading order, namely ei = r/16 and €2 = — r/8 — (n — 1). 
Thus the physical classiflcation of single scalar fleld models is characterized by the 
borderlines €2 = -v^ r = —8{n — 1) and €2 = ei/2 r = —A{n — 1) (short dashed lines 
in Fig.EI). 

This one-to-one correspondence is spoiled by the eg-dependence of n — 1 at the 
second order. Nevertheless, the borderline €2 = (upper full line in Fig. is well 
defined, since higher horizon-flow parameters enter only together with €2 and thus it 
is in principle possible to distinguish toward-exit inflation from hidden-exit inflation in 
the r-tilt plane. 



13 




-0.15 -0.1 -0.05 0.05 0.1 0.15 

n-1 



Figure 2. Classification of models in the r-tilt plane. The vertical long dashed 
line denotes the scale-invariant HZ spectrum. Toward-exit models of inflation have a 
negative tilt (red spectrum) and arc found to the l.h.s. of the upper full line. To the 
r.h.s. of the upper full line we find the hidden-exit models. The (upper) short dashed 
line shows the approximate borderline derived from the leading order expression given 
in the text. The toward-exit models are divided by another thick line which assumes 
vanishing running of the spectral index. Models with decreasing kinetic energy are 
confined between both thick lines. The (lower) short dashed line is the approximate 
leading order result. The tilted long dashed lines show two cases of non- vanishing 
running of the spectral index, namely dn/dlnfc = — 12ef(+4e^), corresponding to 
£3 — -|-2e2(— 262). The upper curve corresponds to the negative sign of the running. 
We also indicate the accuracy of the considered approximations. 

However, the borderline (e2 = ei/2) between the two subclasses of toward-exit 
models, i.e., models with general initial conditions and models with special initial 
conditions, gives rise to a family of lines that depends on the observed running of 
the spectral index, dn/dln/c = — 2eie2 — 6263 (leading order). One possibility is to fix 
63 by assuming that the leading contribution to the running vanishes (lower full line in 
Fig. 12)). Thus it is impossible to distinguish between models with decreasing ("chaotic") 
and increasing ("false vacuum") kinetic energy density on the basis of a r-tilt plot, 
unless the running of the spectral index is known or constrained to be small. 

The above discussion confirms that the horizon-flow functions are more fundamental 
than other quantities that have been used to classify inflation models. However, the 
upper solid line in Fig. |21 is robust in the sense that it holds true for the slow-roll 
approximation and the constant-horizon approximation at any order (as terms e„ with 
n > 2 enter in combination with 62 and are thus zero for €2 = 0). 

Figure |21 shows another remarkable feature of the constant and growing horizon 
approximations: for the shown region in r and n — 1 the amplitudes are accurate to 
0.3%, except for the tiny region enclosed by the arc on the l.h.s. of the figure. This is 
in contrast to the slow-roll approximation at second order, which is less accurate in the 
upper right and lower left corner of the figure. 
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6. Testing the higher-order constant-horizon approximation 

As we already noted, observations show that ei <^ 1. If £2 is of the order of ei or smaUer, 
higher-order corrections (terms that are at least quadratic in the horizon-flow functions) 
to the primordial spectra are irrelevant and the first-order expression for oq, the seed of 
the recursion (|T7jl . is good enough to calculate the coefficients a„ in the power spectrum 
Nevertheless, if €2 is actually much larger than ei, then higher-order corrections in 
€2 could be necessary to match the observational accuracy. In such a case the constant- 
horizon approximation (see Section I4.2.2|) is the simplest and most economic way of 
obtaining high-precision predictions. To see if this is true, let us start by comparing 
Eqs. ()19|) and ()21|). One sees that the third-order expression for the constant-horizon 
approximation is simpler and requires knowledge of less horizon-flow functions than the 
corresponding lower-order slow- roll approximation. Next, we must test how good the 
constant-horizon approximation is. 

To measure the error of the approximations, we define 

AA(l0g ~k) = ^num{l0gk)-A,^^,{\0gk) ^^^^^ ^ ^g^^ 
AnumC^Ogk) 

where A^^^{logk) and ylappr(log A;) stand for the numerical and analytically 
approximated values of the amplitudes at the normalized scale k = k/k^,. 

The constant-horizon approximation applies for many single-field models based on 
phenomenological particle physics: the inverted quadratic model V = Vq — m^(y9^/2 
[5U] . or more generally models with V = Vo[l — {(p/fiY] and p > 2, and also for those 
with V = Vo(l — exp{—(p/ fi)) [H]. For positive fi, all these models belong to the 
class of toward-exit infiation with special initial conditions. There are also examples of 
hybrid infiation models, e.g. those arising from dynamical supersymmetry breaking with 
= Vo[l ± (/i/v?)^], where p is a positive integer |42j. For the positive sign, the model 
belongs to the hidden-exit infiation class, while for the negative sign, the models belong 
to the class of toward-exit infiation with special initial conditions. Yet another model 
belonging to that class is hybrid infiation with a running mass that arises from one- 
loop corrections in supersymmetry-inspired models (see Refs. |^ El 1^ for concrete 
realizations): 

^ - (31) 



V = Vo 



^m2(^)(^2 



Further examples where the constant-horizon approximation applies are the so-called 
natural infiation model : 



V = A' 



1 + cos 




and the model proposed by Wang et al. 



V = A' 



2 , / 

1 arctan (Sy?) 

TT 



(32) 



(33) 
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Figure 3. Test of approximated solutions to the equation of the scalar modes for a 
hybrid model with running mass. Third-order constant-horizon approximation (cha3), 
and first (sral) and second (sra2) order slow-roll approximations are compared. The 
errors shown are those obtained including terms up to 03 in parametrization and 
keeping the highest possible order for the horizon-flow functions in all coefficients. 

which is a designer model that has been used to demonstrate the hmitations of the 
slow-roll approximation. 

We started by testing the natural inflation model given by potential ()H2j) . For 
this model with A = / = 1 and = 1 at the pivot value yj* = 0.01, we obtained 
= 0.81648, ei^ €3^ 10"^ es ^ 0.039529 and €4 ^ 0.02008. With these values, we 
find that the third-order constant-horizon approximation performs slightly better than 
the second-order slow-roll approximation, and that both of them provide a significant 
improvement with respect to the first order slow-roll approximation. The errors are 
confined within the 15% interval in a range exceeding A log k = Q. 

More interesting is a test for a model where some of the higher horizon-flow 
functions are "large". We start from Eq. pip and use the ad hoc choice m'^{(f) = 
M'^ exp {—{p / fi) , with M and fi positive. The resulting potential has a maximum at 
if = and a minimum at = 2fi. Starting near the false vacuum, enough inflation can 
be produced before the minimum is reached. For Vq = 1, fi = 1 and M = 2 at the pivot 
value = 0.01, we find = 0.57729, ei = 0.00003, €2 = 0.29494, €3 = -0.00381, 
64 = 0.14344 and €5 = —0.0082, with the results presented in Figs. IH and El 

As can be observed in Fig. IHl the conclusions drawn for the simple inflation 
near a maximum model are still valid for models where the conditions for constant- 
horizon approximation to apply are met in a weaker fashion than in the case of natural 
inflation, although the range where the errors are confined in the 15% error band is 
smaller. According with the results in Fig. HI increasing the order of the horizon-flow 
functions beyond the quadratic in every included term of parametrization (fTTjl does not 
significantly improve the accuracy of the approximation. 
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Figure 4. Test of the constant-horizon approximation for a hybrid model with running 
mass. The errors shown are those obtained including terms up to a^, in parametrization 
and varying the order of the highest included horizon-flow function. 




Figure 5. Test of the constant- horizon approximation for a hybrid model with running 
mass. Here we investigate the convergence of the Taylor expansion as more and 
more coefficients a„ are added. We keep the highest possible order of the horizon-flow 
functions. 

As shown in Fig. adding terms in parametrization (jllll actually seems to be most 
important for increasing the precision of the prediction, but this conclusion might not 
hold true for other models. Even in the best case, it seems difficult to keep the error 
below the 15% mark for a range broader that A log k = 5. 

We tried to push our approximation to the limits by testing it on the model given by 
Eq. dnSl). For A = 1 and = -0.25, we find = 0.72476, ei = 0.01122, eg = 0.31874, 
63 = 0.19659, €4 = 0.079 and 65 = 0.129. The test confirms the previous results, 
although the range where the error is under 15% is significantly shorter, A log A; ~ 3.5. 

For completeness we note here that similar results were obtained for the tensor 



17 



modes. 

For all tested cases, the constant-horizon approximation at third order performed 
as good or slightly better than the slow-roll approximation at second order, but it has 
the advantage that one needs to make less restrictive assumptions on higher horizon-flow 
functions and the approximation is more efficient in the sense that the recursion seed oq 
only needs ei and €2 as an input, whereas the slow- roll approximation at second order 
needs additionally €3. 
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